On the motion of a heavy rigid body in an ideal 

fluid with circulation 



Alexey V. Borisov, Ivan S. Mamaev 

Institute of Computer Science, Udmurt State University 

1, Universitetskaya str., 426034 Izhevsk, Russia 

Phone/Fax: +7-3412-500295, E-mail: borisov@rcd.ru, mamaev@rcd.ru 

February 8, 2008 

Abstract 

Chaplygin's equations describing the planar motion of a rigid body in 
an unbounded volume of an ideal fluid involved in a circular flow around 
the body are considered. Hamiltonian structures, new integrable cases, and 
partial solutions are revealed, and their stability is examined. The problems 
of non-integrability of the equations of motion because of a chaotic behavior 
of the system are discussed. 

1 



1 Introduction and a review of the known results. 

S. A. Chaplygin jlj considered a general problem regarding the forces and momenta 
that impact an arbitrary rigid body involved in a parallel-plane motion in an un- 
bounded volume of an ideal incompressible fluid. More specific formulations were 
studied earlier by N. E. Zhukovski [221 I2H| , who considered the application of his for- 
mula for the lifting force to the description of heavy body falling in a fluid. However, 
these formulations yielded an unrealistic conclusion that the rotational and trans- 
lational motions are independent, thus showing that additional analysis is required. 
The necessary analysis was made by S. A. Chaplygin. 

Chaplygin made general assumptions that the fluid motion is vortex-free with a 
zero velocity at the infinity, and the fluid circulation around the body is constant. He 
observed that this form of equations holds true if the propeller propulsion does not 
change its direction with respect to the body (aircraft) and the drag equilibrates this 
propulsion at any moment. Although this assumption resulted in the conclusion that 
the aircraft will have a directional stability, which is not true (the aircraft trajectory 
may be winding), some general conclusions that can be drawn from this study allow 
a remarkable mechanical interpretation. They are of both theoretical and practical 
interest even in the context of modern applied aerohydromechanics. 

S.A. Chapygin [3] considered the case of a circulation- free planar flow around 
a body and suggested a remarkable form of equations, which unfortunately he has 
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not studied at all (equations with a similar form describe the falling of a body 
of rotation in a fluid under the same conditions). Originally written in 1890 as 
his student essay, this work was published only in 1933 as part of Chaplygin's 
Collected Works. Independently and at about the same time D. N. Goryachev [S] 
(1893) and V. A. Steklov (1894) obtained the same equations and described 
their simplest qualitative properties. In particular, V. A. Steklov showed that, as 
the body falls, the amplitude of its oscillations with respect to the horizontal axis 
decreases while the frequency of these oscillations increases. He made this re- 
mark in an addendum to his book [201, m which the analysis of the asymptotic 
behavior was made with some mistakes. In essence, V. A. Steklov formulated the 
problem of asymptotic description of the behavior of a heavy rigid body during 
its fall. This problem was solved by V. V. Kozlov j^], who showed that, under al- 
most all initial conditions, the body tends to fall at a uniform acceleration with 
its wider side up and oscillates around the horizontal axis with an increasing fre- 
quency of the order of t and a decreasing amplitude with an order of 1/ y/t. Asymp- 
totic motions with different numbers of half-turns were analyzed numerically in [3] . 
The asymptotic behavior of a body falling without initial impact was studied in 

BB|. 

The effect of an abrupt ascent was described and studied in Under conditions 
of vortex- free flow around the body, it is assumed that, in the initial moment, the 
wider side of the body is horizontal and the body has a horizontal velocity. In the 
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following moments, the body starts moving downward. However, if its apparent 
mass in the lateral direction is sufficiently large, the body will next abruptly move 
upward with its narrower side up and rise higher than its initial elevation. 

More general equations describing a non-planar motion of a heavy rigid body in 
an unbounded volume of an ideal fluid involved in a vortex-free motion and resting 
at the infinity were also obtained in [3| (more compact form of these equations 
is presented in [10J). They generalize the well-known Kirchhoff equations, which 
are known to neglect the force of gravity. As was shown in [TU], if a rigid body 
has three mutually orthogonal planes of symmetry, this body, when falling freely, 
asymptotically tends to take the position when its axis with the maximum added 
mass is vertical. Such body can also rotate around this axis. 

A particular case of existence of the Hessian integral for the general equations [3] 
is described in 0. In this case, the equations can be reduced to a simpler form 
analogous to the planar case. It is worth mentioning that a planar motion of a rigid 
body in a resisting medium is considered in JTj. For various models of hydrodynamic 
forces, some numerical and experimental results concerned with the falling motion of 
a heavy body in a fluid are presented in the papers PQ |2U El H3 HH] • An elementary 
analysis of the motion of a body in a resisting medium was first performed by 
Maxwell in |IE| . 

Now let us return to the problem of a planar circular motion of a rigid body, which 
is the focus of this study. The equations of motion for this problem were suggested 
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by S. A. Chaplygin in [I]. It was studied in [T2| with some specific assumptions (as 
compared with 0]). 

Circulation makes the motion of the rigid body more complicated. In [T2| , cases 
of stationary motion are found and their stability is examined. An integrable case 
and certain classes of particular solutions analogous to Zhukovski's solutions [22H2H] 
are proposed. 

In this paper, we will propose the Hamiltonian form of the general equations, 
and show new cases of integrability and classes of new particular solutions. We will 
also show that in the general case, the equations jl] are not integrable and their 
behavior is chaotic. 

2 General equations of motion. Lagrangian and 
Hamiltonian descriptions. 

We choose a moving coordinate system O1X1X2 which is fixed to the body. The 
position of this system with respect to the fixed frame is characterized by the coor- 
dinates (x,y) of its origin 0\ and the rotation angle ip. Let us assume that (£,77) 
are the Cartesian coordinates of the center of masses in the moving coordinates, I 
is the constant circulation around the body, p is the fluid density; v-i are the 
projections of the linear velocity of the center 0\ onto the moving axes. Let us 
assume also that Qi, Q 2 are the projections of the external forces onto the moving 
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axes, oj is the angular velocity of the body, and Mi is the momentum of the external 
forces relative to the center of masses. 

Explicit evaluation of the forces and momenta associated with the circulation 
flow around a body was used in j3j to derive the following equations, which are 
similar to the Kirchhoff equations: 

divi = cLiLuv 2 - Xv 2 — (uj — Qi, 
a 2 v 2 = -a 2 0JVi + Xvx + - Q2, 
boo = (ai - a 2 )v 1 v 2 + (>i - X v 2 - M 1 + £Q 2 - tjQi, 

(1) 

(f> = W, 

x = (a, v) = V\ cos <p — v 2 sin ip, 
y = (7, v) = vi sin ip + v 2 cos <p, 
where a. = (cos^ — siny?), 7 = (sia<p, cosy?), v = (vi,v 2 ). 

In equations (JTJ) a\,a 2 ,b are added masses and momentum, A = Ip, ( = Ip\i\, 
X = Ip\i- Parameters (, Xi proportional to the circulation J, are associated with 
the asymmetry of the body, and their evaluation is described in 0] . In the general 
case with circulation and asymmetry of the body ( / 0, x 7^ 0, these parameters 
cannot be eliminated by any choice of a coordinate system fixed to the body. 

Remark. Equations (0) were proposed by S.A.Chaplygin in 1926 and somewhat 
later (and presumably independently) were derived by Lamb and Glowert in 1929 
\14j - Unlike Chaplygin ^j, the behavior of solutions of these equations was not 
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analyzed in these works. 

It is worth mentioning that equations similar to (JTJ describe the planar motion 
of a rigid body in an ideal fluid with a uniform vorticity. The motion of a circular 
cylinder in a fluid with a uniformly distributed vorticity (for different boundary con- 
ditions) was analyzed by Prudeman and Taylor In this case, the equations of 
motion are reduced to a simplest linear system. The general (even planar) motion 
of an arbitrary rigid body is very complicated and has not been analyzed yet. 

If all forces are potential, we have 

dU dU . 

Ql = cos Lf + — sin cp, 
ox ay 

Q 2 = — q — sin ip H - — cos (p, {*) 
ox oy 

dU 

Mi - CQ 2 + rjQi = 

where the function U = U(x, y, ip) is the potential. In this case, equations (JHJ) have 
the integral of energy 

H = T+U = -{aivt + a 2 vt + ko 2 ) + U, (3) 

where the kinetic energy T of the ( "body + fluid" ) system is written in the diagonal 
form after translations and rotation of the moving coordinate system. 

It can be shown that the forces caused by circulation are generalized potential 
forces and the equations of motion (in the case of potential external forces) can be 
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written in the Lagrangian form 



d (9L\ dL dL dL 

dt \ dv\ ) dv 2 dx 1 dy 

d ( dL\ dL dL dL 

dt \ dv2 ) dvi dx 2 dy ^ 2 

dt \du J V<2 dvi 1 dv 2 dip 



(4) 



<P = w, x= (a,v), y = (j,v), 
where the Lagrangian function has the form 

L = T — U — (x sin (p — y cos (p) (^ v i + X^J ~ ( x cos V 9 + V srn V 9 ) (^ V2 + • (^) 

Remark. The Lagrangian function is chosen in a symmetric calibration; recall 
that the summands that are linear with respect to velocities are determined to a total 
differential of an arbitrary function. 

Equations are Poincare equations on the group of motions of the plane 
E(2) 2\; using the Legendre transformation we can represent them in the Hamilto- 
nian form (Poincare-Chetaev equations) with a Hamiltonian containing terms linear 
in the momenta. It was found , that the equations of motion for generalized po- 
tential systems can be conveniently represented in the slightly modified variables: 

dT dT „, <9T 7 

pi = — = a x vi, p 2 = — = a 2 v 2 , M = ^—buj. 
av\ dv 2 duo 

Using the Legendre transformation for new variables, we find the Hamiltonian: 



The obtained equations of motion have the form similar to that of the Poincare- 
Chetaev equations: 

dH dH dH dH dH 

Pi = - OLl ~ -7T- 7l - ^ C" 



' dM dx dy dp 2 dM 

dH dH dH dH dH 
^ 2 ^ dM dx 012 dy ^ 2 dpi ^ dM ' 

dH dH dH dH dH 

M = Pi a P 2 ^ ~a~ + C"» ' 

dp2 dpi dip dpi dp2 

dH . ( dH\ . ( dH\ 
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v= m> x= \ a >-dj)> v= K t >-dj) 

and the Poisson bracket of these variables contains "additional circulation terms" 
(hygroscopic terms): 

{Ml, pi} = -p 2 + C, {M,p 2 } = Pi-Xi {PhPz} = -A, 

{M,<p} = -l,{M,x} = {M,y} = 0, (8) 

{p h x} = -a 1} {p 2 ,x} = -a 2 , {pi,y} = -71, {P2,y} = -72- 
The rank of this Poisson structure is equal to 6; therefore, the system (jZJ can be 
reduced to a canonical system with three degrees of freedom. 



3 Motion in the gravity field. 

In this case, the potential energy of the system can be written in the form U = 
+ £ sin (p + 77 cos <p), and its Hamiltonian (jjj) has the form 

1 (p\ p\ M 2 \ 
H=- (— + — + -T- )+ /i{y + £ an <p + v cos <P)- (9) 



2 V ai 02 & 



The system (J7j). Q admits an autonomous and a non-autonomous integral, which 
correspond to the projections of the momentum on the fixed axes Ox and Oy (see 
figure 1): 

Px = Pi cos ip — P2 sin if + Xy + ( sin if — x cos ip = c\ , 
P y — P2sinyj + p 2 cos(p — Xx — (cosip — X sin if = /it + c 2 , (^) 
Ci = const, c 2 = const, 
which are not commutative {p x ,Py} = X. 

To perform a reduction by one degree of freedom we express y from p x = C\ and 
substitute it into the Hamiltonian 

f/ = 1 H ) + -r (-Pi COS f+p 2 Sin <p) + 



cos y? 



This reduced Hamiltonian depends only on the variables p\, p 2 , M and </?, the 
Poisson brackets of which, according to (jHJ), form a closed subalgebra with a rank 
of 4. Thus, we have a reduced system with two degrees of freedom, which can be 
written in the canonical form (see below). However, we will not use the canonical 
form, but algebrise the reduced system to even greater extent with the help of 
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dependent variables 71 = simp, 72 = cos<^ : 

i =1(- + - + t)-i«'-* , '> + ** 

9M ^ (12) 

p-s,3—j -/i(r,J 7 ), 

. _dH 

1 ~ m J7, 



where p = (pi,p 2 ), 7 = (71,72), J 



1 



s = (a,C), r = (£,77). 



-1 

The system has two obvious integrals — energy H and geometric integral 7J + 
7§ = 1. Since the system is Hamiltonian, for being integrable it requires one more 
additional integral (although this system can be integrated by the Euler-Jacobi 
method because the equations (JT2j) conserve the standard invariant measure). 

Note that the system (JT2j) has an important mechanical meaning, its form is as 
simple as that of the Euler-Poisson equations, and analogous formulations of prob- 
lems are quite meaningful for this system. The first aspect relates to the integrability 
of the system ([12)1 ■ 

Canonical variables. Using the algorithm described in [2], it is easy to construct 
canonical variables analogous to the Anduaye variables in the rigid body dynamics. In 
conventional denotations, we have 



M = L, pi - x = y/2X(G- L)costp, p 2 - C = y/2X(G - L) sinZ, 

M + fr-* )a + <*-< )a = G , y = B , 
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where I, L and g, G are two pairs of canonical conjugate variables ({Z,L} = {g,G} = 1). 
These variables can be useful for construction of action-angle variables and for application 
of methods of the Hamiltonian theory of perturbations. 



4 Integrab le cases. 



The following integrable cases are known. 



1. n — (S. A. Chaplygin, 1926 [4J — the absence of the gravity field. 



Additional quadratic integral has the form 



(13) 



and the system (|12|) can be reduced, similar to the classical Euler-Poinsot case, to 
a system of three equations for variables M, p±, and p 2 - Chaplygin in jl] noted that 
explicit integration involves some very complicated quadrature, which is expressed 
in elliptic functions under the condition \ = £ = 0. 

2. C = x — 0, a\ = a 2 = 1 (V. V. Kozlov, 1993 [12]) — a case of dynamical 
symmetry. This case is analogous to the Lagrange case, although the additional 
integral is quadratic with respect to momenta 



F 



1M 2 
2~b~ 



+ /z(r,7). 



(14) 



Let us consider a new case of integrability. 
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3. Let us assume that x = ~A?7, C = A£ (i- e - s = — AJr). Additional quadratic 
integral has the form 

F = iM 2 -A6(M + (p,Jr))+^(r, 7 ). (15) 

Except the cases mentioned above, this system doesn't have other (general) cases 
with an additional integral linear or quadratic with respect to phase variables. This 
simple statement can be proved by direct enumeration of variants with the use of the 
method of undetermined coefficients. The question of the existence of higher-order 
integrals is still an open question. 

Let us consider a linear invariant relation analogous to the Hess case for the 
Euler-Poisson equations. 

4. Let us suppose that ( = x — 0, r\ = & ( a 2 1 ~ a :f 1 )- m this case, 

2 \ai a 2 b J A V 
F = M±sJb (a 2 x -a^ 1 ) 

and at the level F = we have F = 0. 

Remark. A^oie that in the case with zero circulation (A = 0), the reduction to 
the system (|T2*|) is impossible since j = oo. Nevertheless, integrals p x — 7Ti ; p y — 
7r 2 + fit, 7Ti, 7r 2 = const allow us to eliminate p\, p 2 from the equations (rather than 
from the Hamiltonian) and obtain, similar to Chaplygin J3]/ ; the non- autonomous 
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"pendulum" -type equation for <p: 



p. t , i. , n at ( , _i _ lw . _ . xcosv 9 C sirL 9^ \ 

p = — ;— (a 2 — a i ) sin Zip + — ( a 2 — % J ( Ci cos 2(p + c% sin 2(p) 1 

2b b \ a\ a>2 J 

— cosy? — rjcosp) + — — ^ 2 - ((c 2 — c 2 ) sin2<p — 2c 1 c 2 cos2p) 

p . 7Ti / xsinp Ccos(/? 
— — (4 cos <p — 77 sm p) + — I h 



b b \ a\ a 2 

7r 2 / xcosy? C srn ¥> 



b \ d\ 0,2 

(17) 

which in particular cases was examined in [3J, 



5 Chaplygin's case. Bifurcation analysis. 

Let us consider the system (JUJ) at p = 0, p = x = 0. As was noted above, equations 
for M, pi, P2 can be separated and have the form 

Pi = p 2 M - A—, p 2 = — + A—, M — pip 2 . (18) 

a 2 a 2 \a 2 a x J 

The common level of the first integrals 

H= \(rt + Pl + ^_\ =h = const; F = M + = c = const (19) 

2 \ai a 2 b J 2A 

is represented by closed curves formed at the intersection of an ellipsoid and an 
elliptic paraboloid; these curves are analogous to centroid lines in the classical Euler- 
Poinsot problem (figure 2). It is easy to find and describe the particular solutions 
of (J18)) which are analogous to permanent rotations in the latter problem. 
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For the sake of definiteness, let us assume that A > and ai > a 2 . 

I. pi — P2 — 0, M — c — const, h = ^c 2 . In this case, the body rotates 
uniformly (<p = | = const) around the origin of the moving coordinate system 0\ 
(see figure 1), though this origin is fixed (and in the general case does not coincide 
with the center of mass). This solution is unstable when — < c < — and stable 
otherwise (the stability on the boundaries requires special analysis). 

II, III. There exists a pair of analogous solutions of the form 

p. = 0, M = — b = const, p) = 2A ( c - — ) , (20) 

where in one case, % — 1, j — 2, while in the other case, % — 2, j — 1. Each of these 
solutions exists under condition c > — , and the constants of the integrals for them 
are related by 

h = c --b^. (21) 
aj aj 

In the case of such motion, the solid body uniformly rotates around the origin of 
a fixed coordinate system, and the principal axes of the body always pass through 
this fixed center. The body is directed toward the center of rotation by its wider side 
in one motion and narrower side in the other. The solution p 2 = 0, h = ^-c — \^*zb 
is always stable while the solution p 1 — 0, h — -^c — \ ^*zb is always unstable. Thus, 
stable motions are such that the body is directed toward the fixed point by its 
narrower side. 

The bifurcation diagram is shown in figure 3. The straight lines represent the 
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particular solutions II and III; they are tangent to the parabola corresponding to 
solution I in the points where admissible values of the integral c begin. The domain 
of possible motions on the bifurcation diagram is hatched. 

The bifurcation diagram presents three different intervals of the values of energy h 

for which the patterns of the trajectories corresponding to different values of c are 
qualitatively similar (figure 2). 

Let us consider explicit formulas for bi-asymptotic motions for solution I, which 
is unstable at — < c < — . These solutions are homoclinic (see figure 2 b) and have 
the form 

2Bd 1 d 2 



M 



cosh 2r — d\ + d\ 



2 ' 



/8A5 djdlcoshr 8\B rff^sinhr , . 

lh ~" d\ C osh2r -dl + d 2 2 ,P2 ~ V d\ cosh2r-d? + ^' 1 " ? 



c . , „ / 1 + d? - d\ , \ v / M 2 " 

V= £(*-*o) -2arctan I — — tanhr I , r= — 

where df = |, d\ = |, 6i = 2(c - A&a^ 1 ), 6 2 = 2(A6a^ - c), £ = &i + 62, and 
different signs correspond to the two different separatrices. 
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6 Perturbation of Chaplygin's case. Splitting of 
separat rices. 

Let us consider a "perturbed" Hamiltonian , H = H + with /i regarded as a 
small parameter. Here, H is the Hamiltonian of the integrable Chaplygin problem, 
Hi = U is the perturbation function. One of the dynamic effects that prevent the 
existence of an additional analytical integral for the perturbed system is the splitting 
of separatrices, which are coupled in the case of the unperturbed system. To the 
first order of the theory of perturbations at small //, the splitting of separatrices is 
determined by the value of the Poincare-Mel'nikov integral 



calculated along asymptotic solutions of the unperturbed problem. Here F is an 
integral of the unperturbed system. 

If the integral (as a function of the parameter on the separatrix) has a simple 
zero, then separatrices split and transversally intersect. In the case of two degrees 
of freedom this makes the perturbed problem non- integrable [T5] . 

We take the integral (JED F = M 2 + ^1 as p . 

On homoclinic asymptotic solutions (|22|). the integral ()23j) can be explicitly eval- 




(23) 
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uated using residues: 



J{go) = Jq + Jr, 



Jq = -2tt 



\/2/i6 sinh (| - a) S 



sm g Q , 



(24) 



Jr = -4vr 




(£, cos g - n sin g ) 



i? sinh vrS* 



where S = 2c(& 1 5 2 )~ 1 , cos a = d\, sin a = c? 2 - 

The integral (J2*%j) has a simple zero, except for the special case of B = 0. This 
allows us to conclude that the perturbed system (jHJ) is non-integrable, except for 
the case of a 1 = a 2 . Note that it is under this condition that the integrals (JHj), (|T5|) 
were obtained. This condition is necessary but not sufficient for the integrability. 

The chaotic behavior of the system (jHJ), associated with the splitting of separa- 
trices and non-integrability, is illustrated in figure 4, which represents a Poincare 
section of the phase flow of the system (fT2"|) on the energy level of H = const. The 
figure shows the behavior of the split separatrices and the stochastic layer that forms 
near them. The figure also presents the projection of the section at the energy level 
of H = h = const, which is determined by the relationship ip = it = const, on the 
plane (pi,^)- As before, the parameter values are taken as follows: A = 1, b = 2, 
a\ = 2<i2 = 1, x = C = an d for the energy level at these parameter values we 
assume that h = 2. 

Let us discuss the issue of the possible fall of a body in the presence of circulation. 
Recall that the body moving in the absence of circulation asymptotically tends to 
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fall with its wider side down [§]• 

// A ^ 0, the moving body will stay in a finite-width band parallel to the x axis. 
Indeed, according to (fTU|). we have 

x = -jt + A -1 (Oi - x) sin<f + (p 2 - C) cos<p - c 2 ) , 
V = -A -1 ((pi - x) cos cp - (p 2 - C) sin if - d) ( 25 ) 
Ci, c 2 = const, 

We assume that pi, p 2 are limited and rewrite (flip. The right-hand part of this 
equality is clearly a limited function at A 7^ 0, therefore p±, pi are also limited. 

According to (|25|). the body moves in the horizontal direction with a mean 
velocity of — ^. This result was obtained in [T3] for x — C — 0- 
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Figure 1: 



a) h = 1/2 6) h = 3 c) h=7 

Figure 2: Trajectories of the system in Chaplygin's case in the space of variables p%, 
p<2, M. The parameter values are: A = 1, b — 2, a\ — 1, 02 = 31 X — C — 0. 
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ai > a 2 




Figure 3: Bifurcation diagram for Chaplygin's case at a\ > a 2 . The domain of 
physically admissible values of integrals c, h is shown in gray. Unstable permanent 
rotations are shown by dashed lines. Parameter values are: A = 1, 6 = 2, a\ — 
1, a 2 = \. 
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Figure 4: Transition to chaos in the Chaplygin system with an increase in //. The 
first three figures show the trajectories and one pair of separatrices, while, the sep- 
aratrices are not shown in the latter two figures because the respective fixed point 
becomes stable at given parameter values. (The domain of variables where the mo- 
tion is impossible at the given parameter values is shown in gray. The scale of the 
last figure differs from that of the others.) 
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